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Abstract 

In this paper, we investigate some sufficient conditions for the breakdown of local smooth 
solutions to the three dimensional nonlinear nonlocal dissipative system modeling electro¬ 
hydrodynamics. This model is a strongly coupled system by the well-known incompressible 
Navier-Stokes equations and the classical Poisson-Nernst-Planck equations. We show that the 
maximum of the vorticity field alone controls the breakdown of smooth solutions, which reveals 
that the velocity field plays a more dominant role than the density functions of charged particles 
in the blow-up theory of the system. Moreover, some Prodi-Serrin type blow-up criteria are 
also established. 
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1 Introduction 


In this paper, we study a well-established mathematical model for electro-diffusion, which is gov¬ 
erned by transport, and Lorentz force coupling between the Navier-Stokes equations of an incom¬ 
pressible fluid and the transported Poisson-Nernst-Planck equations of a binary, diffuse charge 
system (cf. [29]). The three dimensional Cauchy problem reads as 


dtu + {u ■ V)u — fj,Au + Vn = eAiE'ViE', 

X G t > 

V • u = 0, 

X G t > 

dtv + {u ■ V)u = V ■ (UiVt) — vivV"^), 

a; G t> 

dtw + {u ■ \7)w = V • {D 2 yw + r' 2 tcV'I'), 

X G K^, t > 

= V — w, 

X G K^, t > 

{u,v,w)\t=o = {uo,Vo,'Wo), 

X G K^, 


*E-mail addresses: jihzhao@163.com (J. Zhao), 124204602@qq.com (M. Bai). 
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where u = {u^,u'^,v?‘) is the velocity vector field, 11 is the scalar pressure, is the electrostatic 
potential, v and w are the densities of binary diffuse negative and positive charges (e.g., ions), 
respectively, is the kinematic viscosity, e is the dielectric constant of the fluid, known as the 
Debye length, related to vacuum permittivity, the relative permittivity and characteristic charge 
density. Di, £> 2 , vi, V 2 are the diffusion and mobility coefficients of the charget0. We refer the 
reader to see m for the detailed mathematical description and physical background of this fluid- 
dynamical model. For the sake of simplicity of presentation, we shall assume that all physical 
parameters to be one in the system (HH). 

The system HU was introduced by Rubinstein in |29j . which is intended to account for the 
electrical, fluid-mechanical and bio-chemical phenomena simultaneously occurring in complex bio- 
hydrid system like bio-chips, bio-reactors or micro-fluidic chambers in Lab-On-Chip technology, 
see nnni Elisa for specific applications of the system HU and nananu for the computational 
simulations. Generally speaking, the self-consistent charge transport is described by the Poisson- 
Nernst-Planck equations, and the fluid motion is governed by the incompressible Navier-Stokes 
equations with forcing terms. 

To our knowledge, mathematical analysis of the system HU was initiated by Jerome m, 
where the author established a local well-posedness theory of the system HU based on the Kato’s 
semigroup framework. Concerning existence results of weak solutions, we refer the reader to see 
na ini i3Qi EH, where the system HU adapted with various boundary conditions are considered 
(e.g., Neumann, no-flux, mixed boundary conditions), while considering existence results of smooth 
solutions, we refer the reader to see [UEUEH], where local existence with any initial data and global 
existence with small initial data in various functional spaces (e.g., Lebesgue, Besov spaces) are 
established. Since the Navier-Stokes equations is a subsystem of HU, one can not expect better 
results than for the Navier-Stokes equations. Hence, in the case of three dimensional space, whether 
the regularity and uniqueness of global weak solutions or global existence of smooth solutions are 
still the challenge open problems. Some regularity and uniqueness issues have been studied by the 
references mimiiiiH] even for more general system for the electro-kinetic fluid model. 

System HU is closely connected to classical models in fluid mechanics and in the theory 
research of electrolytes. When the charge transport is dropped, i.e., r, = re = 'k = 0, one gets the 
well-known three dimensional incompressible Navier-Stokes equations: 

{ dtu -I- (u ■ V)u — /rAu -I- VH = 0, x G K^, t > 0, 

V-u = 0, xGM.^, t>0, (1.2) 

u\t=o=uo, a; e 

This equations has drawn considerable attention of researchers for many years, and many fruitful 
results can be found from the literature, e.g., lUllinilllllUlMlESlES] and the references therein. 
In the pioneering work |24j , Leray proved global existence of weak solutions with any initial data in 
in dimensions two and three, moreover, the global existence of smooth solutions in dimensions 
two and in dimensions three under smallness condition imposed on the initial data were also 

^Di = D 2 = , where Tq is the ambient temperature, k is the Boltzmann constant, and e is the 

charge mobility. 
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addressed. It is well-known that in dimensions three, without smallness condition imposed on 
initial data, whether the corresponding local smooth solutions can be extended to the global one is 
an outstanding open problem (alternative open problem is the regularity and uniqueness of global 
weak solutions). The well-known Prodi-Serrin criterion (cf. [351131]) shows that any weak solution 
u for the Navier-Stokes equations satisfying 

u £ L^(0,T;LP(]R^)) with --f - = 1 and 3 < p < oo, (1.3) 

q p 


then u is regular on (0,T] x Beirao da Veiga |4] extended the condition (11.311 to the following 
condition: 


Vu G L'^(0, T; L^(]R^)) with —|-- = 2 and ^ < p < oo. 

q p 2 


(1.4) 


In 1984, Beale, Kato and Majda in their celebrated work |3] showed that if the smooth solution u 
blows up at the time t = T, then 


[ ||a;(-,t)||Lcx, dt = oo, 

Jo 


(1.5) 


where w = V x u is the vorticity. Subsequently, Kozono and Taniuchi |23j and Kozono, Ogawa 
and Taniuchi [22], respectively, refined the criterion (lESD to the following two criteria: 


riH- 

Jo 


\bmo dt = oo and 




B<L 


dt = < 


( 1 . 6 ) 


where BMO is the space of Bounded Mean Oscillation and B^ oo(®^) is the homogeneous Besov 
spaces. Recently, Dong and Zhang [5] improved the criteria (fra to the Navier-Stokes equations 
(II.2p that if 

f \\VhU^i-,t)\\BO dt<oo, (1.7) 

Jo 

then the solution u can be extended beyond the time T, where V/j { 81 , 82 ), and {v},vJ) 

is the horizontal components of the velocity field u. 

Based on the basic energy inequalities, by using the standard fixed point argument, it is easy 
to show local well-posedness of the system (HID for sufficiently smooth initial data and hence we 
have the following result. 


Theorem 1.1 Let uq £ with V • uq — 0, vq, wq £ T^(IR^) n and vq, wq > 0. Then 

there exists a time T* > 0 such that the system (HD admits an unique solution {u,v,w) satisfying 

w £ C'([0, T], i72(R3)) n L°°(0, T; n L‘^{ 0 , T; H^{R^)) 

for any 0 < T < T*. Moreover, v,w >0 a.e. in [0,T*) x 
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Inspired by [3], the first purpose of this paper is to establish an analog of Beale-Kato-Majda’s 
criterion for singularities of local smooth solutions to the system (HD. The result shows that the 
maximum norm of the vorticity alone controls the breakdown of smooth solutions, which reveals 
an important fact that the velocity field u plays a more dominant role than the charge density 
functions v and w in the blow-up theory of smooth solutions to the system p.ll) . We follow 
closely the arguments used in [3] to accomplish the proof in the following three steps: obtaining 
estimates for the vorticity w and v,w, obtaining higher energy estimates for the solution {u,v,w), 
and applying the crucial logarithmic Sobolev inequality. 

The main results of this paper are as follows: 

Theorem 1.2 Let uq £ with V ■ uq = 0, vq^wq G L^(R^) fl and vq^wq > 0. Let 

Ti, > Q he the maximum existence time such that the system (HD admits an unique smooth solution 
(u, u, w) satisfying (11.811 for any 0 < T < T*. IfT^< oo, then 

\\uj{-,t)\\L’^ dt = oo. (1.9) 

In particular, 



limsup ||a;(-,f)||L“ = oo. 

t/'T, 

Remark 1.1 As a byproduct of the proof of Theorem ll.il we obtain that under the assumptions 
of Theorem 11.21 if for some 0 < T < oo such that 

||a;(-,t)||Loo dt < 00, 

then the solution (u, v, w) can be extended past the time t = T. 

Motivated by the blow-up criteria (HD-(HD for the Navier-Stokes equations (HD, the second 
purpose of this paper is to generalize these blow-up criteria to the system HD- 



Theorem 1.3 Let uq G with V ■ uq = 0, vo,wo G L^fU.^) n and vo,wo > 0. 

Assume that T* > 0 is the maximal existence time such that the system (HD admits an unique 
solution (u, V, w) satisfying (HD for any 0 < T < T*. Then T, < -|-oo if and only if one of the 
following conditions holds: 


\\u{-,t)\\‘fpdt = oo with 


H— =1 and 3 < p < oo, 
P 


( 1 . 10 ) 


\\'Vu{-,t)\\fpdt = oo with —I—=2 and 


P 


3 
2 

n 


< p < oo. (1-11) 

and Vo, wq > 0. 


Theorem 1.4 Let uq G with V • uq = 0, uo,wo G L^{ 

Assume that T* > 0 is the maximal existence time such that the system (HD admits an unique 
solution (u,v,w) satisfying (HD for any 0 < T <T„. Then T* < -|-oo if and only if the following 
condition holds: 


rT, 


WVhU^i- 


\lo 


2p 

P’-r 


2 3 3 

dt = oo with —I— —2 and — < p < oo, (1-12) 


q p 


where Vh = {di,d 2 ) and 
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Remark 1.2 It is clear that when p = oo, the condition (11.121) becomes jJ'* || po dt = oo, 

thus Theorem 11.41 can be regarded as an extension of the blow up criteria (fra-lfTTl). and is even 
new for the Navier-Stokes equations (O. 

Before ending this section, we recall the definition of the homogeneous Besov spaces, for more 
details, see [1] and [34]. Let cj) : ^ [0,1] be a smooth cut-off function which equals one on the 

ball |) '^= G : |^| < |} and equals zero outside the ball |). Write 

p - (j){2^) and 

Then we define for all j G Z, 

A,/ p{2-W)f. 

We have the following formal decomposition: 

f = J2A,f for fGS'{R^)/V{M.^), (1.13) 

iez 

where is the set of polynomials (see [1]). Let sGR, l<p,r<oo and / G 5'(R^), the space 

of temperate distributions, the norms in homogeneous Besov spaces are defined as follows: 

dg I for 1 < r < oo, 

\supjgz2J"||Aj/||LP forr = oo. 

Then the homogeneous Besov space i?p_^(R^) is defined by 

• For s < - (or s = - if r = 1), we define 

{/ G 5'(R3) : < CX)}. 

• If A: G N and ^ + k < s < |-|-A:-|-1 (ors = |-|-fc-|-lifr = l), then i3p_^(R^) is defined as 
the subset of distributions / G 5'(R^) such that f G 5'(R^) whenever |/3| = k. 

Another classical results that will be used here are the so-called Bernstein’s inequalities (cf. 
my- there exists a constant C such that for every / G <S'(R^), any nonnegative integer k and any 
couple of real numbers (p, q) with 1 < p < g < oo, we have 

supp/C {|e| < 2^} ^ sup WIWl. < (1.14) 

\(y.\—k 

supp/c{|e|«2^} ^ C-i-'=2^'=||/|Up< sup ||a“/|Up<Ci+'=2^'=||/|Up, (1.15) 

\ot\—k 

where / represents the Fourier transform of /. 

The rest of this paper is organized as follows: We shall present the proof of Theorem 11.21 in 
the Section 2, while the proof of Theorem 11.31 will be presented in Section 3 and Theorem O in 
Section 4. Throughout the paper, we denote by C and Ci {i = 0,1,---) the harmless positive 
constants, which may depend on the initial data and its value may change from line to line, the 
special dependence will be pointed out explicitly in the text if necessary. 
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2 Proof of Theorem 11.2 


We prove Theorem 11.21 by contradiction. It suffices to show that if the condition (EJl) does not 
hold, i.e., there exists a constant K > 0 such that 


then one can prove that 



||u;(-,<)||Lood< < K, 


( 2 . 1 ) 


sup (||■u(•,^)||//3 + < Co, (2.2) 

0<t<T, 

where Cq is a constant depending only on ||mo||h 3, ||('yoj^fo)||Linff2j 21* and K. This implies that 
the solution (u, v, w) can be extended beyond the time T*, which contradicts to the maximality of 
T,. 

Let us first recall the following two important inequalities that play an important roles in the 
proof of Theorem 11.21 for its proof, see [20] and [3] , respectively. 

Lemma 2.1 ([20]) If s > 0 and 1 < p < oo, then 

WA^ifg) - fA^gh. < C{\\AfU^\\A^-^gU. + (2-3) 

where A = (—A )2 , and C is a constant depending only on s and p. 

Lemma 2.2 ([3]) For all u € 7L^(M^), w = V X It, there exists a constant C such that 

||Vu||l- < C(1 + ||a;||L 2 + ||w||loo ln(e + ||n||ff3)). (2.4) 

Next we establish the bounds for the vorticity w, v and w under the assumptions of Theorem 
o and the condition (EU. 


Lemma 2.3 Assume thatuo G iL^(K^) withV-uo = 0, {vo,wo) S L^(]R^)nL^(K^) andvo,wo > 0. 
Let (u, V, w) be the corresponding local smooth solution to the system HH) on [0, T) for some 
0 < T < 00. If 


||w(-,t)||L~ dt = K < OO, 

then we have 

sup (\H;t)\\h + M-Mh) + 2 r (l|Vn(.,t)||i. + ||Vn;(.,t)||i.)dt < Co 
0 <t<T ^ ' Jo ^ ^ 

and 

sup ||w(-,t)||i2 + / ||Va;(-,t)||| 2 dt < Co, 

0<i<T Jo 



(2.5) 


( 2 . 6 ) 


(2.7) 


where Cq is a constant depending only on ||ito||i?i; ||(uo,u;o)||LinL2, T and K. 
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Proof. We first recall from m that if Vo and wq are non-negative, then v and w are also non¬ 
negative. Next, in order to derive bounds of v and w, we introduce two symmetries: ^ ‘^= v + w 
and 77 = V — w. Then we infer from (ED that C and 77 satisfy the following system: 


'9tC+(u-V)C = AC-V-(77V«'), 
dtT] + {u ■ V)77 = A 77 — V • (CV'I'), 
A'h = 77 , 

^Co = ^'o + Wo, rjo = vo -wq. 


Multiplying the first equation of (12.81) by C and integrating over one has 
^4llClli2 + ||VC||i.+ / CV77-VvI/dx+/ C^Ad7dx = 0, 
where we have used the fact 

f {u-V)C(^dx = ^ [ u-V((^)dx = 0 
Jr^ 2 Jg3 


/R 3 

due to V • 77 = 0. The third equation of 
1 d 


R3 

yields that 


2 dt 


IICII 


L2 


+ l|VC|li 2 + / (:Vrj-V^dx+ [ C,rfdx = Q. 
JR3 7r3 


It can be done analogously for 77 , we get 

^4:lhlli3 + + / 77VC-V^'dx-f / Cv^dx = 0. 

^ ® dR3 dR3 

Adding (12.91) and ( 12 . 101 ) together, after integration by parts, we obtain 


l_d 

2 dt 


(llCIlia -I- ||77|||2^ -b ||VC||i 2 -I- IIV 77 III 2 -b C,rfdx — 0. 


2 

L2, 


(2.9) 


( 2 . 10 ) 


( 2 . 11 ) 


Integrating (12.111) over [0,t] for all 0 < t < T implies that 

llC(0lli3 + h(0lli2 +2^‘ (l|VC(r)||i 2 + ||V77(T)||i2)dr 

+ 2 / [ Cix,T)r]^[x, T)dxdT = \\Co\\l 2 + 

Jo dR3 

which gives us to 

\\v{t)\\l2 + \\w{t)\\l2 + 2 (||Vw(t)||| 2 -b ||V77;(T)|||2)dT 

+ 2 / / iv{x, t) + w{x, t)){v{x, t) - w{x, T)fdxdT = \\vo \\\2 + ||wo||i 2 . 
Jo Jr^ 

Since v and w are non-negative, we infer from (12.121) that for all 0 < t < T, 

\\v{t)\\l2 + ||w(t)||i2 +2 f (||Vv(T)||i2 -b ||V77;(T)|||2)dT < ||7;o||i2 -b ||wo||i2. 


( 2 . 12 ) 
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This proves (1^ . 

To derive bound of u, multiplying the first equation of (HID by u, after integration by parts, 
by the fifth equation of (EH), it can be easily seen that 

+ l|VM||i2 = (2.13) 

On the other hand, we multiply, respectively, the third and fourth equations of EH by 'P to yield 
that 


J + V ■ + (u ■ dx = 0, 

J (^dtw'^ — V • (rcV'I')'!' — Aw^E' + {u ■ V)w'i’^dx = 0. 

Subtracting (12.151) from (12.141) leads to 

J — w)4' + V • ((w + w)V'I')4' — A{v — w)'^/ + (m • V)(v — dx = 0, 
which, after integration by part and using the fifth equation of EH again, gives us to 
i^||V4'|||2 + [ (t! + w)|V4fpcia: + f |A4'pcia:+ f {v — w)u ■V'i’dx = 0. 

2 dt J^3 J^3 Jr3 

Adding (12.131) and (j2.16|) together, we find that 

■{Mh + IIVVE-Ili^) + W^uWh + \\Anh +Jj^ + = 0. 


l_d 
2 dt' 


Observing that v and w are non-negative, thus we have 


1 d 

2 dt 


u\\h 


IV^I 


h) + IIVmII 


2 

L2 


|A4/||i2 < 0. 


Integrating (12.171) in the time interval [0, t] for any 0 < t <T, one obtains that 


<t)ll 


L2 


l|V4'(t)||i2+2 


(2.14) 

(2.15) 


(2.16) 


(2.17) 


r)||i2 + ||AvI/(r)||i2 dr < ||uo|li. + llVvhoHi^, (2.18) 


where 4'o is determined by vq and wq through the equation A^'o = vg — wg. By using the Sobolev 
embedding relation ^ L®(]R^), the Young inequality and the interpolation inequality, we 

obtain that 


||V4'o||i2 < ||4 'o||l6||A4'o|1^6 < ||4'o||l 6||(uo,wo)||^6 

< C'||V4'o||l2|I(^o,'Wo)||^ 6 < -||V4'o|li2 + C||(uo,w;o)||% 

< ^l|V4'o|li2 -f C'||(^;o,M;o)|IEl||(^'o,wo)|||2• 

Hence, for all 0 < t < r, we obtain the following bound of u and V'l': 

Mmh + l|V'P(t)||i 2 + 2 f (||Vu(T)||i 2 + ||A4^(r)||i2)dr < WugWh + Ci, 

^0 


(2.19) 












where Ci is a constant depending only on ll(t'O) w'o)||LinL 2 . 

Finally, we establish the desired estimate for w. Taking Vx on the first equation of CH) , we 
see that 

dtuj + {u ■ V)a; — Aw = (w • V)u + V x (A'I'V'I'). (2.20) 

Multiplying (12.201) by w, and integrating over one has 


2 dr 




/R3 


(A'kV'I') • (V X oj)dx 


def j 
= 


( 2 . 21 ) 


By the Biot-Savart law, we have the following two equalities: 

du 

— = Rj{R X w) for j = 1,2,3, 


and 


V X w = V X (V X u) = V(V ■ u) — Au = —Au, 

where R = (i?i,i? 2 ,.Ra), Rj = (j = 1)2,3) are the Riesz transforms. Since the Riesz 

operators are bounded in L^(R^), we have 

||V'u||l 2 < C'||w||i 2 and ||Au||i 2 < C'||Vw||i 2 . 


This implies that 


h < 



{lu ■ \7)u ■ Lodx 


<||c^||L<»||Vu|U2||w|U2<C||a;|U)»||a;||i.. 


Moreover, by using the fifth equation of (fTTl) . (1^ and (12.191) . we bound I 2 as 


/2 <c||v^'|U4||A4-|U4||Am|U2 

<\\\Au\\l.+C\\{v,w)\\l 4Vnl^ 

< l||Vw||i.+C||(u,w)||i||(Vr;,Vrc)||{.||Vvk|||,||Avk||{, 

< l||Vw||i.+C||(u,w)||i.||(Vr;,Vu;)||i||VvI-|||, 

< l||Vw||i.+C||(u,w)||i.||(Vr;,Vrc)||i.+C'||(^;,w)||i.||VvI/||i. 
<i||Vw||i.+C(||(V^,Vu;)||i.+l), 

where we have used the following interpolation inequality: 

||/||L4<||/||}.||V/||f.. 

Plugging (|2.22l) and (12.231) into (12.211) . we obtain 

^Il<^lli2 + l|Va;|||2 < C^IIwIIloo + ||(V?;, Vw)|||2 + 1^ (^||w|||2 + 1^. 


( 2 . 22 ) 


(2.23) 


(2.24) 
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Applying the Gronwall’s inequality, we conclude that for all 0 < t < T, 




L2 


||Va;(T)|||2dT < Caexplc^ (||a;(T)||L=o + ||(Vv(t), Vu;(t))||^2 + l)dT|, (2.25) 


where C 2 = ||wo|||2 + 1- By (12.51) . (12.61) and (12.251) . we get (12.71) . The proof of Lemma is 
complete. □ 

Now, let us derive higher order derivative bounds of the solution (u,v,w). Taking VA on the 
first equation of (HH), multiplying the resulting with V An and integrating over , observing that 
the pressure If can be eliminated by the condition V ■ n = 0, one obtains 

~\\VAu\\l 2 + \\A'^u\\l 2 = - [ VA((u • V)it) • VAudx + f VA(A^'V«') • VAudx 

2 dt Jjja 7 r3 


Ih+Ih. 


(2.26) 


Since V • u = 0, it follows from Lemma 12.II that 
Ih = - 


/ [VA((u • V)u) — {u ■ V)VAn] • VAudx 

< C'||VA((n • V)u) - (u ■ V)VAu||L2||VAn||L2 

< C'||Vn||Loo||VAn||i2. 


(2.27) 


For 7/2, by using the fifth equation of (irD) . (EH), EH) and (j2.19|) . one has 
II 2 = — [ A{{v — w)V'^/) ■ A^udx 

JR? 

\\AMl^+C\\A{{v-w)Vml^ 

\\A^u\\l 2 + c(||(Au - An;)V«'||i2 + 2||(Vn - Vw)VH\\l 2 + ||(n - n;)VA^'||2,) 
||A2n||i2+c(||(Au,An;)||i3||Vd/||ie + ||(n,n;)||i3||(Vn,Vn;)||ie) 

\\A^u\\l 2 + C^IKAu, Au;)||i2||(VAn, VAu;)||i2 + ||(Vu, Vw)||l 2 ||(An, An;)||^2 


1 

< - 

- 2 
1 

< - 

- 2 
1 

< - 

- 2 
1 

< - 

- 2 
1 

< - 

- 2 


|A2n||i2 + i||(VAn,VAn;)||i2+c(l+||(Vn,Vn;)lli2)||(An,An;)||i2, (2.28) 


where we have used the interpolation inequality: 


Il3 < 


lL2 


iiv/iii2, 


and the following Hardy-Littlewood-Sobolev inequality with p = 2: 

l|V(-A)“Vll ^<C\\f\\LP for any 1 < p < 3. 

L 3-p 

Putting (12.271) and (12.281) together, we deduce that 
|||VAn||i. + IIA^nlli^ < i||(VAn, VAn;)||i3 
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+ c(l + II VulUoo + IKVz;, VHIIi.) (|| VAwlli. + IKAz;, Aw)\\l.). (2.29) 

Taking A to the third equation of (ll.ip . then multiplying Av, after integration by parts, we see 
that 

-^||Aw ||^2 + ||VAn ||^2 = — [ A{{u ■ \7)v)Avdx — f AV ■ {v'V'^)Avdx 
2 dt 7r3 jr3 

+ J//2. (2.30) 

Applying Lemma [2731 the terms IIIi and III 2 can be estimated as follows: 

///i = / V((it • V)v) ■ VAvdx 

dR3 

< i||VAn||i2 + C(||(Vw • V)r;||i2 + ||(u • V)Vr;||i2) 

< i||VAn||i 2 +C(||V«||i„„||Vr;||i 2 + ||«||ie||AnHia) 

< i||VAn||i2 + C(||Vn||i2(l + ||VAw||i2) + ||Vu||i2||A^;|U.||VA^|U3) 

< -||VAw ||^2 + (7(1 + ||Vn||^ 2 ) (1 + ||VAm ||^2 + ||An||^ 2 ) 

and 

III 2 = I A(z)V'I') • S/Avdx 

dR3 

< i||VA^;||i 2 + C(||AWVI.||2, + ||VW2d/||i2 + |kVAvI/)||i.) 

< i||VA^;||i2 + C(||A^||i3||Vd/||ia + ||V^;||ie||V2vl,|||3 + ||^||i31|V(^; - «;)||ia) 

< i||VA^||2, +C(1 + IKVt;, Vzz;)||i2)(l + ||(Az;, Ara)||i2), 

which give us to 

j^WAvWl. + ^||VAn||i2 < C(1 + IKV^;, Vz/;)||i2)(l + ||VAu||i2 + ||(Az;, A«;)||i2). (2.31) 

It can be done analogous for w, thus we get 

|||Au;||i2 + ^||VAu;||i2 < C(l + ||(Vz;, Vu;)||i2) (l + ||VAu||i2 + ||(Au, Au;)||i2). (2.32) 

In the final step, we complete the proof of Theorem 11.21 bv applying the logarithmical Sobolev 
inequality (12.21) and the Gronwall’s inequality. Set 

Y{t) = e+Wumj,. + \\{vit),wm\l^- 


Then by (12.61) . (I2.19L (12.291) . (12.311) and (12.321) . we see that 

^ < ^(l + ||Vu|U<» + ||(Vu, Vu;)||i2)F(t). 


(2.33) 
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On the other hand, using and the logarithmic Sobolev inequality in Lemma 12.21 we get 


||Vu||loo < 0(1 + IIwIIloo ln(e + (2.34) 

Combining (I2.33|) . (12.341) with the fact Inm(t) > 1 yield that 

^lnr(t) <C'(l + ||a;|U=. + ||(Vu,Vu;)||i.)lnr(t). (2.35) 

By applying the Gronwall’s inequality, we obtain 

y(t) < exp I InF(O) exp ((7(1 + t) + C J ||a;(r)||L=o(ir) | (2.36) 

for any 0 < t < T*, where F(0) = e + Huoll^a + ll(V0,^i'0)||ff2. Hence, by (I2.36|) . we obtain 

sup (||u(-,t)||^/3 + ||(u(-,t),-u;(-,t))||/f 2 ) < Co, 

0<t<T. 

where Co is a constant depending only on ||uo||//3j IK'i'O!'u^o)||Lini? 2 , 7* and K. The proof of 
Theorem o is complete. 


3 Proof of Theorem 11.3 


In this section we prove Theorem ll.3l Based on the Beale-Kato-Majda’s criterion (11.91) in Theorem 
11.21 we need only to show that if one of the conditions ( 11 . 101 ) and ( 11 . 111 ) does not hold, then there 
exists a constant C such that 

\\u{-,t)\\H3dt < C. (3.1) 

By the Sobolev embedding L°°{R^), (13.11) gives us to the fact \\uj{-,t)\\Lc^dt < C, 

which completes the proof of Theorem 11.31 bv the Beale-Kato-Majda’s criterion (11.91) in Theorem 

O 



Lemma 3.1 Assume thatuo € withV-uo = 0, {vo,wo) G L^(]R^)nL^(K^) andvo,wo > 0. 

Let {u, V, w) be the corresponding local smooth solution to the system dD on [0,T) for some 
0 < T < oo. If there exists a constant K > 0 such that one of the following conditions holds: 


then we have 


2 3 

/ WuLjCWlpdt < K with — I — = 1 and 3 < p < oo, 

Jo <1 P 

/■^ 2 3 3 

/ \\Vui-,t)\\%pdt < K with — I —=2 and — < p < oo, 

Jo g P 2 


sup \\yu{-,t)\\l 2 + [ \\Au{-,t)\\l 2 dt < Cs, 

0<t<T Jo 


(3.2) 


(3.3) 


(3.4) 


where Co is a constant depending only on ||uo||i?i; IK'i'O: ii'o)||LinL3, T and K. 
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Proof. Multiplying the first equation of (ED by Au, after integration by parts, we deduce that 


- —||Vm |||2 + ||Au |||2 = — [ u-Vu-Audx+ [ ■ Audx‘^= Ji + J 2 - (3-5) 

2 dt 7r3 Jr3 

Under the assumption (13.21) . using the Holder inequality, Ji can be estimated as 

Ji <C||u||lp||Vm|1^^||Au|U 2 <C||M||Lp||VM||^;"||AM||^t" 

<\\\Au\\l.+C\\u\\U^u\\l., (3.6) 

where q = while under the assumption (13.31) . Ji can be done as 
^' r 

Ji = / dku^diudkudx < C\\Vu\\lp\\'^u\\'^ ^ 

< C\\Vuh4^u\f~J\\Au\\l < ^\\Au\\l,+C\\Vu\\UVu\\l., (3.7) 

where q = 2 p -3 ■ •^ 2 ) by using the fifth equation of (11.11) . (12.61) and (I2.19L we get 

J2<c\\v^\\l4^^\\l4M\l- 

<^\\Au\\i.+c\\{v,w)\\i4ynh 

< i||Au||i2+C||(^,u;)|||,||(Vu,Vu;)|||,||VvI/||J,||A^-|||, 
<i||Au||i2+C||(^,u;)||i2||(Vu,Vu;)|||,||VvI/||J, 
<i||Au||i2+C||(^,u;)||i2||(Vu,Vu;)||i2+C||(^,u;)||i2||VvI/||i2 

< i||Au||i2+C(||(Vu,Vzc)||i2 + l). (3.8) 

Plugging p.6p and p.8p into (13.51) . we have 

+ I|Am||| 2 < C'||u|||p||Vu|||2 + CdKVu, Vw;)|||2 + 1) 

<Ci\\u\\% + ||(V^, Vu;)||i 2 + 1) (e+ ||Vu||i 2 ) . (3.9) 

Applying the Gronwall’s inequality yields that 

sup ||VM(f)||i 2 < (e+||Vuo||i 2 )exp j / (||M(r)|||p + ||(Vu(r), Vw;(t ))|||2 + l)dr| , 

0<t<T [JO J 

which combining p.6l) . (12.191) and (j3.9l) give us to (|3.4I) . 

On the other hand, plugging (13.7p and (13.8p into (|3.5I) . we have 

|||Vzr||i 2 + ||AM||i 2 < C (llVwIlL + ||(Vu, V«;)||i 2 + l) (e + ||Vu||i 2 ) . (3.10) 
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Applying the Gronwall’s inequality again yields that 


sup ||Vu(t )||^2 < (e+ ||Vuo||^ 2 )expi f (||Vu(t)|||p + ||(Vu(r), Vw(t))||^ 2 + l)dT i , 
0<t<T [JO J 

which still leads to (j3.4p by (I2.6p , (j2.19|) and (I3.10|l . The proof of Lemma 13.11 is complete. □ 

Now we prove (EH). Multiplying Av to the third equation of (II.IL integrating over we 

obtain 

-^||V '(;|||2 + ||Au |||,2 = — [ {u-V)vAvdx— f V • (uV4')Auda;. (3-11) 

2 at 7r3 7r3 

Applying Lemma IXTl the two terms in the right hand side of (13.101) can be estimated as follows: 
- f (u ■ V)vAvdx < i||Au||i 2 + C||(u • V)v\\l 2 

jR3 O 

A gl|Au||^2 + C'llulllellVullla 

<i||Au||i2+C'||Vu||i2||Vu|U3||Au|U3 
- + C'||Vu|||2, 


- / V • {vV'i>)Avdx < i||Au||i 2 + C'(||Vz;V4'||i2 + ||uA4'||i2) 


A gl|Aw |||2 + C'dlVulllallV^'llle + II(u,re) 111,4) 

<i||Au||i2+C(||(u,u;)||i2||Vu|U3||Au|U3 + ||(u,u;)||j2||(Au,Au;)||J.) 

<i||(Au,AHIli3+C(l + ||(Vu,VHIli3). 

Notice that same result also holds for w, thus we obtain 

^(||Vu||i2 + ||Vu;||i2) + ^(||Au||i2 + ||Au;||i2) < C(1 + ||(Vu, Vu;)||i2). (3.12) 


Taking A on the first equation of dlH, multiplying the resulting with Au and integrating over 
by the condition V • u = 0, we see that 

- —||Au |||2 + ||VAu |||2 = — [ A{{u ■'V)u) ■ Audx + ( A(A'I'V4') • Audx. (3.13) 

2 dt 7r3 JR3 

Applying Lemma [XT] again, the right hand side of (13.131) can be bounded as follows: 


I A{{u ■ V)m) • Audx < -||VAm |||2 

R3 8 

<^l|VAu||i2 

<^l|VAu||i2 


+ C||(Vu-V)u||i2 + ||(u.V)Vu||i2 
+ C||Vu||l4 + ||u||ia||Au||i3 
+ C||Vu||i.||VAu||f, + ||Vu||i2||Au|U2||VAu|U. 
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<-||VA«||i. + C(||A«||i. + l), 

/ A(A'I'V'I') • Awda; = / V((z; — u>)V'I') ■ VAwda; 

Jr3 Jr3 

< IwVAuWh + c(||(Vz; - Vw)Vnh + ll(^ - 

<i||VAu||i.+c(||(Vt>,Vzi;)||i3||VvI/||ie + ||(z;,zz;)||i4) 
<i||VAu||i. + i||(A^,Au;)||i.+c(l + ||(V^,Vzz;)||i3). 

Hence, we have 

j^W^uWl. + \NAu\\l. < i||(An, A«;)||i. +C'(l + ||(V^, Vu;)||i. + ||An||i.). (3.14) 

Finally, we conclude from (13.121) and (13.141) that 

|(||Au||i. + ||(Vn, V«;)||i.) + (||VAzr||i. + ||(Az;, An;)||i3) 

<C(l + ||(Vn,Vn;)||i3 + ||Au||i.), 

which yields directly (13.11) by the Gronwall’s inequality. We complete the proof of Theorem 11.31 


4 Proof of Theorem 11.41 

In this section we prove Theorem ll.4l Based on the proofs of Theorems ll.2l and ll.31 it is sufficiently 
to show that if there exists a positive constant K such that 


2 3 3 

I • „ dt < K with — I — =2 and - < p < oo, 
vM <i P 2 


then dSl still holds, which gives us the desired estimate m- Actually, we shall establish 
bounds of the vorticity w = V x u, which leads to the desired bounds of Vu by the boundedness 
property of Riesz operators in L^. Taking Vx on the first equation of m, we see that 


dtLo + {u ■ V)a; — Aw = (w • V)u + V x (A'f'VT). 
Multiplying (14.111 by w, and integrating over R^, one has 

;||i 2 + ||Vw||i 2 = / {uj-V)u-ujdx- [ (AW^-) • (V x a;)da 

Jr3 Jr3 


(4.1) 


Id 

2 dt' 




(4.2) 


Notice that 


Ki = / tUidiU^ojjdx = / tUidiU^ojjdx + / uiidiU^ui^dx 


ij=l 


/R3 


ij=l ' 


i=l ■ 
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Thus, by using the facts 



ijJ^d^u^wzdx. 


/R3 


= -{div} + d2V?), W 3 = diu^ - 92^^, 


we can easily see that 

Ki < f \VhU^\\i^\'^dx. 

JK3 

Applying the Littlewood-Paley dyadic decomposition (11.131) for i-e., 


N 

^ AjVhU^ + ^ AjVhU^ + ^ AjVhU^, 

j<-N i=-N j>N 


we can split Ki into 


(4.3) 



Now we estimate the terms Kij [j = 1,2,3) one by one. For ATn, by using the Holder inequality, 
the Bernstein inequality (I1.14|) and the fact ||VM||i 2 < C'||a;||i 2 , one has 

Kii = [ \Aj'^hU^\\u}\‘^dx 

i<-Af 

<c ^ \\A,y,u^U^Ml. 

3<-N 

j<-N 

<c( ^ ^ ||A,V,n'^||i.)^|a;||i. 

j<-N j<-N 

<C2-^^\\yu\\Ml. 

<C2-i^Mh. (4.5) 

For Ki 2 , let p' be a conjugate index of p, i.e., ^ + A = 1 . it follows from the Holder inequality 
that 

N 

Ki 2= y] / \AjVhU^\\ujfdx 

j=-N 

N 

<c y \\A,S/^u'^U4oj\\l,^, 

j=-N 

N 2p-3 N 2 

<c( £ 1^)^( £ \\A,Vnu4t)^’’Ml.p 

j=-N j=-N 
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(4.6) 


P’^ 

<l\\Vu;\\h+CN\\VHU%, ^ Mh, 

P'T 

where q = 2 p- 3 ; have used the following Gagliardo-Nirenberg inequality: 

ll/ll^^ <||/|ll;^||V/||f. for p>|. 

For Ki 3 , notice that 


IIV^ulU. = ||Au|U 2 <C||Va;|U2. 

Therefore, based on the above inequality, applying the Holder inequality, the Bernstein inequalities 
(I1.14|) and (11.151) again, K 13 can be estimated as 



< C ^ \\AjVhU^\\L3\\uj\\L^\\w\\Le 

3>N 

< C ^ 25||A,V;,u'*|U2||a;|U2||Va;|U2 

j>N 

<C[Y. 2”0'( E 2^^||A,V,u^||i.)"||.;|U.||Va;|U2 

j>N j>N 

<C2-^\\uj\\l2\\V\\\l2\\Vuj\\l2 

<C2-^\\ujU2\\S/u:\\1.. (4.7) 

Plugging estimates (IrSt^diTl) into dSH) gives us to 

K, < hvujWl. + C2-i^Ml, + CN\\S/hu%, Ml. + C2-^||c.|U2||Va;||i.. (4.8) 

8 p.¥ 

Notice that K 2 has been treated as I 2 in (j2.23|) . thus inserting estimates (14.81) and (12.231) into (14.21) . 
we obtain 

+-MujWI. < C2 Wu^Wl. + CN\\y^ \M\l'^+C2 2 ||a;||^ 2 ||Vw |||2 

P'W 

+ C(||(Vu,Vu;)||i2+l). (4.9) 

Now if we choose N sufficiently large such that 

C2 2 ||w||^2 < 

i.e., 

ln(||a ;||?2 + e) + 21nC 

N > - - -h 2, 

In 2 
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then the inequality (14.91) yields that 

+ l|Vw |||2 < ^(l + ||V?iu''||^o ^ + ||(V?;, Vw)||^ 2 ) (11^11^2 + e) In (||a ;||^2 + e) . (4.10) 

Setting Z{t)‘^= e + ||w||^ 2 - It follows from (14.101) that 

|lnZ(t) < C{1 + \\yHU%, ^ + ||(V^, Vu;)||i 2 ) InZ(f). 

Integrating on the time interval [0,t], then taking maximum norm with respect to t for any 0 < 
t < T,, we conclude that 

rT, 

sup Z(t) < Z(0) exp I exp {C / (l + || + ||(Vw(t), Vw(T))||| 2 )dr}|. 

0<t<T, '' Jo p,^ 

Namely, 

sup ||w( 0 lli 2 < (e+||wo||i 2 ) 

0<t<T* 

X explexpICT,+ (:; / (||V,iu'*(t)|||o + ||(Vu(t), Vri;(T))|||2)(ir}|, 

Jo p,^ 

which gives us to m G L°°(0, T*; D L^(0, T,; i?^(]R^)). We complete the proof of Theorem 

[El- 

Acknowledgments. J. Zhao is partially supported by the National Natural Science Foundation of 
China (11501453), the Fundamental Research Funds for the Central Universities (2014YB031) and 
the Fundamental Research Project of Natural Science in Shaanxi Province-Young Talent Project 
(2015JQ1004). M. Bai is partially supported by the Foundation for Distinguished Young Talents 
in Higher Education of Guangdong (2014KQNCX223). 


References 

[1] H. Bahouri, J.-Y. Chemin, R. Danchin, Fourier Analysis and Nonlinear Partial Differential 
Equations, Grundlehren der mathematischen Wissenschaften, vol. 343. Springer, Berlin, 2011. 

[2] P. B. Balbuena, Y. Wang, Lithium-ion Batteries, Solid-electrolyte Interphase, Imperial Gollege 
Press, 2004. 

[3] J. T. Beale, T. Kato, A. Majda, Remarks on breakdown of smooth solutions for the 3D Euler 
equations. Comm. Math. Phys. 94 (1984) 61-66. 

[4] H. Beirao da Veiga, A new regularity class for the Navier-Stokes equations in R", Chin. Ann. 
Math., Ser. B 16 (1995) 407-412. 

[5] J. Bourgain, N. Pavlovic, Ill-posedness of the Navier-Stokes equations in a critical space in 
3D, J. Funct. Anal. 255 (2008) 2233-2247. 


18 




[6] L. Caffarelli, R. Kohn, L. Nirenberg, Partial regularity of suitable weak solutions of the Navier- 
Stokes equations, Comm. Pure Appl. Math. 35 (1982) 771-831. 

[7] C. Deng, J. Zhao, S. Cui, Well-posedness for the Navier-Stokes-Nernst-Planck-Poisson system 
in Triebel-Lizorkin space and Besov space with negative indices, J. Math. Anal. Appl. 377 
(2011) 392-405. 

[8] B. Dong, Z. Zhang, The BKM criterion for the 3D Navier-Stokes equations via two velocity 
components. Nonlinear Analysis: Real World Applications 11 (2010) 2415-2421. 

[9] H. Fujita, T. Kato, On the Navier-Stokes initial value problem. Arch. Rational Mech. Anal. 
16 (1964) 260-315. 

[10] E. T. Enikov, B. J. Nelson, Electrotransport and deformation model of ion exchange membrane 
based actuators. Smart Structures and Materials 3978 (2000) 129-139. 

[11] E. T. Enikov, G. S. Seo, Analysis of water and proton fluxes in ion-exchange polymer-metal 
composite (IPMC) actuators subjected to large external potentials. Sensors and Actuators 
122 (2005) 264-272. 

[12] J. Ean, H. Gao, Uniqueness of weak solutions to a nonlinear hyperbolic system in electrohy¬ 
drodynamics, Nonlinear Anal. 70 (2009) 2382-2386. 

[13] J. Fan, F. Li, G. Nakamura, Regularity criteria for a mathematical model for the deformation 
of electrolyte droplets. Applied Mathematics Letters 26 (2013) 494-499. 

[14] J. Fan, G. Nakamura, Y. Zhou, On the Gauchy problem for a model of electro-kinetic fluid, 
Applied Mathematics Letters 25 (2012) 33-37. 

[15] J.W. Jerome, Analytical approaches to charge transport in a moving medium, Tran. Theo. 
Stat. Phys. 31 (2002) 333-366. 

[16] J.W. Jerome, The steady boundary value problem for charged incompressible fluids: 
PNP/Navier-Stokes systems. Nonlinear Anal. 74 (2011) 7486-7498. 

[17] J.W. Jerome, R. Sacco, Global weak solutions for an incompressible charged fluid with multi¬ 
scale couplings: Initial-boundary-value problem. Nonlinear Anal. 71 (2009) 2487-2497. 

[18] L. Jin, J. Fan, The vanishing viscosity limit for a 3D model of electro-kinetic fluid in a bouned 
domain. Applied Mathematics Letters 26 (2013) 154-157. 

[19] T. Kato, Strong solutions of the Navier-Stokes equations in M'" with applications to weak 
solutions. Math. Z. 187 (1984) 471-480. 

[20] T. Kato, G. Ponce, Commutator estimates and the Euler and Navier-Stokes equations. Comm. 
Pure. Appl. Math. 41 (1988) 891-907. 

[21] H. Koch, D. Tataru, Well-posedness for the Navier-Stokes equations, Adv. Math. 157 (2001) 
22-35. 


19 



[22] H. Kozono, T. Ogawa, Y. Taniuchi, The critical Sobolev inequalities in Besov spaces and 
regularity criterion to some semi-linear evolution equations, Math. Z. 242 (2002) 251-278. 

[23] H. Kozono, Y. Taniuchi, Bilinear estimates in BMO and the Navier-Stokes equations. Math. 
Z. 235 (2000) 173-194. 

[24] J. Leray, Sur le mouvement d’un liquide visqueux emplissant I’espace, Acta Math. 63 (1934) 
193-248. 

[25] M. Longaretti, B. Chini, J.W. Jerome, R. Sacco, Electrochemical modeling and characteriza¬ 
tion of voltage operated channels in nano-bio-electronics. Sensor Letters 6 (2008) 49-56. 

[26] M. Longaretti, B. Chini, J.W. Jerome, R. Sacco, Computational modeling and simulation of 
complex systems in bio-electronics, J. Computational Electronics 7 (2008) 10-13. 

[27] M. Longaretti, G. Marino, B. Chini, J.W. Jerome, R. Sacco, Computational models in nano- 
bio-electronics: simulation of ionic transport in voltage operated channels, J. Nanoscience and 
Nanotechnology 8 (2008) 3686-3694. 

[28] G. Prodi, Un teorema di unicita per le equazioni di Navier-Stokes, Ann. Math. Pura Appl. 48 
(1959) 173-182. 

[29] 1. Rubinstein, Electro-Diffusion of Ions, SIAM Studies in Applied Mathematics, SIAM, 
Philadelphia, 1990. 

[30] R. J. Ryham, Existence, uniqueness, regularity and long-term behavior for dissipative systems 
modeling electrohydrodynamics, arXiv:0910.4973vl. 

[31] M. Schmuck, Analysis of the Navier-Stokes-Nernst-Planck-Poisson system. Math. Models 
Methods Appl. Sci. 19(6) (2009) 993-1015. 

[32] J. Serrin, On the interior regularity of weak solutions of the Navier-Stokes equations, Arch. 
Rational Mech. Anal. 9 (1962) 187-195. 

[33] M. Shahinpoor, K. J. Kim, Ionic polymer-metal composites: 111. Modeling and simulation as 
biomimetic sensors, actuators, transducers, and artificial muscles, Smart Mater. Struct. 13 
(2004)1362-1388. 

[34] H. Triebel, Theory of Function Spaces, Monogr. Math. vol. 78, Birkhauser Verlag, Basel, 1983. 

[35] B. Wang, 111-posedness for the Navier-Stokes equations incritical Besov spaces Adv. 

Math. 268 (2015) 350-372. 

[36] T. Yoneda, 111-posedness of the 3D Navier-Stokes equations in a generalized Besov space near 
BMO-\ J. Funct. Anal. 258(10) (2010) 3376-3387. 

[37] J. Zhao, C. Deng, S. Cui, Global well-posedness of a dissipative system arising in electrohy¬ 
drodynamics in negative-order Besov spaces, J. Math. Physics 51 (2010) 093101. 


20 


[38] J. Zhao, T. Zhang, Q. Liu, Global well-posedness for the dissipative system modeling electro¬ 
hydrodynamics with large vertical velocity component in critical Besov space. Discrete Contin. 
Dyn. Syst. Ser. A 35(1) (2015) 555-582. 


21 



